Double universality of a quantum phase transition in spinor condensates: 
the Kibble-Zurek mechanism and a conservation law 
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We consider a phase transition from antiferromagnetic to phase separated ground state in a spin-1 
Bose-Einstein condensate of ultracold atoms. We demonstrate the occurrence of two scaling laws, for 
the number of spin fluctuations just after the phase transition, and for the number of spin domains in 
the final, stable configuration. Only the first scaling can be explained by the standard Kibble-Zurek 
mechanism. We explain the occurrence of two scaling laws by a model including post-selection of 
spin domains due to the conservation of condensate magnetization. 
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The idea of a noncquilibrium phase transition has been 
attracting great attention in many branches of physics. 
The number of physical models where it has been consid- 
ered or observed is steadily growing, now including not 
only the dynamics of the early Universe [l [and superfluid 
Helium Q , but also superconductors J] , cold atomic 
gases [Ej], and other systems (fjj. The most notable out- 
come of such a phase transition is the possible creation 
of defects, such as monopoles, strings, vortices or soli- 
tons [H 0, 0] after crossing the critical point at a finite 
rate. The Kibble-Zurek mechanism (KZM) is a theory 
that allows to predict the density of created defects from 
the knowledge of the correlation length £ at the instant 
when the system goes out of equilibrium Q • The result- 
ing scaling displays universal behaviour and is dependent 
only on the critical exponents of the system v and z. 

Quantum phase transition, in contrast to a classical 
(thermodynamic) one, occurs when varying a physical 
parameter leads to a change of the nature of the ground 
state 0] . Recently, a few theoretical works demonstrated 
that the KZM can be successfully applied to describe 
quantu m p hase transitions in several models see 
Ref. hj for reviews. Among these, Bose-Einstein 
condensates of ultracold atoms offer realistic models of 
highly controllable and tunable systems 12, 13 1. In fl3j | . 



the miscibility-immiscibility phase transition leading to 
the formation of stable, stationary domains was proposed 
as an ideal candidate to observe KZM scaling in a quan- 
tum phase transition. 

In this Letter, we investigate the critical scaling of the 
number of defects created during the transition to the 
phase separated state of an antiferromagnetic spin-1 con- 



densate [15Hl7j . The introduction of a weak magnetic 



field can lead in these systems to the transition from an 
antiferromagnetic ground state to a state where domains 
of atoms with different spin projections separate [lj|. In 



contrast to the transition considered in [131 ] . our scheme 
does not require the use of microwave coupling field or 
Feshbach resonances, which makes the experiment sim- 
pler and more stable against inelastic losses. 



By employing numerical simulations within the trun- 
cated Wigner approximation, we make a quite unex- 
pected observation. While the number of spin fluctua- 
tions just after the phase transition closely follows the 
predictions of the KZM, the number of spin domains in 
the final, stabilized state is given by a scaling law with 
a different exponent. We explain this double universal- 
ity and predict the value of the second exponent using a 
model of system dynamics including later stages of evo- 
lution, no longer described by the standard KZM. We 
show that when the nonlinear processes set in, an ef- 
fective post-selection of spin domains takes place, after 
which only a part of them can survive due to the conser- 
vation of the condensate magnetization. 

We consider a dilute antiferromagnetic spin-1 BEC in 
a homogeneous magnetic field pointing along the z axis. 
We start with the Hamiltonian H = Hq + Ha, where the 
symmetric (spin independent) part is 



H a = 



E 

J=-,0,+ 



Co 



2m 



where the subscripts j = — ,0,4- denote sublevels with 
magnetic quantum numbers along the magnetic field axis 
m.f = — 1,0, +1, m is the atomic mass, n = = 
J^Wj^j is the total atom density and V(x) is the ex- 
ternal potential. Here we restricted the model to one 
dimension, with the other degrees of freedom confined 
by a strong transverse potential with frequency lo± . The 
spin-dependent part can be written as 



H A = 



J dx ^£> i + |:F 2 :j 



where Ej arc the Zccman energy levels and the spin 
density is F = 'F x ip,i/)^ 'Fyip,ij;^ F z ip) where F x , y , z 
are the spin-1 matrices and ip = (ip+jipo,^-)- The 
spin-independent and spin-dependent interaction coef- 
ficients are given by Co = 2hu±(2a2 + ao)/3 > and 
C2 = 2hu)_i_(a,2 — &o)/3 > 0, where as is the s-wave scat- 
tering length for colliding atoms with total spin S. The 
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FIG. 1. Ground state phase diagram of an antiferromagnetic 
condensate for magnetization M — N/2. We increase B lin- 
early during the time tq to drive the system through one or 
two phase transitions into a phase separated state. 



total number of atoms N = J ndx and magnetization 
M = J (n+ — n_) dx are conserved. 

The linear part of the Zecman effect induces a ho- 
mogeneous rotation of the spin vector around the di- 
rection of the magnetic field. Since the Hamiltonian is 
invariant with respect to such spin rotations, we con- 
sider only the effects of the quadratic Zeeman shift [l6| . 
For sufficiently weak magnetic field we can approximate 
it by a positive energy shift of the to/ = ±1 sublevels 
5 = (E+ + E- - 2E )/2 w B 2 A, where B is the magnetic 
field strength and A = (gi + gj) 2 /^/16-Ehfs, gi an d 
gj are the gyromagnctic ratios of electron and nucleus, 
fiB is the Bohr magneton, -Ehfs is the hyperfine energy 
splitting at zero magnetic field |16j . 

The ground state phase diagram, shown in Fig. [lja), 
contains three phases divided by two critical points at 
B i = B Q M/V2N and B 2 = B Q /V2, where B Q = 
\J cinj A and n is the total density. The ground state 
can be (i) antiferromagnetic (2C) with ip = (ip + ,0,ip_) 
for B < Bi, (ii) phase separated into two domains of 
the 2C and ip = (0,^0,0) type (p ) for B e {B 1 ,B 2 ), 
or (Hi) phase separated into two domains of the po 
and ip = (V>+,0,0) type (p+) for B > B 2 0. Note 
that the phase diagram is different in the special cases 
M = 0, ±iV, which we disregard in the current paper. 

For simplicity, we consider a system in the ring-shaped 
quasi-lD geometry with periodic boundary conditions at 
±L/2 and V(x) = 0. The magnetic field is initially 
switched off, and the atoms are prepared in the antifer- 
romagnetic (2C) ground state with magnetization fixed 
to M = N/2 (without loss of generality). To investigate 
KZM we increase B linearly during the time tq to drive 
the system through one or two phase transitions into a 
phase separated state. Due to the finite quench time, the 
system ends up in a state with multiple spin domains. 

The Kibblc-Zurek theory is a powerful tool that allows 
to predict the density of topological defects resulting from 
a nonequilibrium phase transition without solving the full 
dynamical equations. The concept relies on the fact that 
the system does not follow the ground state exactly in 
the vicinity of the critical point due to the divergence of 
the relaxation time. The dynamics of the system cease 
to be adiabatic at t ~ — t (here we choose t = in the 
first critical point), when the relaxation time becomes 



7"rcl 



where e(t) = B — B\ ~ t/rq is the distance of the system 
from the critical point. At this moment, the fluctuations 
approximately freeze, until the relaxation time becomes 
short enough again. After crossing the critical point, dis- 
tant parts of the system choose to break the symmetry in 
different ways, which leads to the appearance of multiple 
defects in the form of domain walls between domains of 
2C and po phases. The average number of domains is 
related to the correlation length £ at the freeze out time 



t 



ZV J (1 + Zf) 



Na 



(2) 



where z and v are the critical exponents determined by 
the scaling of the relaxation time t, 



rcl 



\~ zv and exci- 
tation spectrum ui ~ |fc| z , with z = 1 in the superfluid. 

We test the above prediction in numerical simulations 
within the truncated Wigncr approximation, with ./V = 
10 6 and parameters close to that of previous experiments 



in d Na [15(. Additionally, we consider a "cleaner" case 
with transition through the first critical point only and 
N = 20 x 10 6 in order to minimize merging of domains. 
Typical results of a single run, which can be interpreted 
as a single realization, are shown in Fig.[2j We can clearly 
see the process of domain formation after the first phase 
transition at t\ . However, there is always some number of 
spin fluctuations that disappear instead of evolving into 
full domains, see Fig. Eff) . 

The above dynamics have a striking effect on the 
number of defects that are created in the system. In 
Fig. [3] we show the average number of defects in the 
function of the quench time tq for the "cleaner" case 
of large number of atoms, as in Fig. EJe). The crit- 
ical exponents, calculated from the Bogoliubov excita- 
tion spectrum of the relevant gapped mode with Ab = 
c^nsj (5/c2n — l) 2 — (1 — M 2 /N 2 ) are identical as in the 
case of a ferromagnetic [l2[ or two component [l3| con- 
densate, z = 1 and v = 1/2. According to the formula 

— 1 /3 

([2]), we could expect the scaling Nd ~ Tq . However, 
the number of domains in the final, stabilized state scales 

2/3 

approximately as N^ ~ t q m a wide range of tq, as 
depicted with the dashed line. Nevertheless, if we count 
the number of spin fluctuations [18j just after the first 
phase transition, we do recover the Nf ~ Tq depen- 
dence (inset). 

We explain this puzzling appearance of two different 
scaling exponents with a model that includes several 
phases of the domain formation process, see Fig. |4j The 
2C state becomes unstable shortly after crossing the first 
critical point at B = B\. Initially, the system follows 
the standard Kibble-Zurek scenario, as the spin fluctua- 
tions start to grow cxponantially at the freeze out time 
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FIG. 2. Spin domain formation dynamics in a ring-shaped ID geometry with ring length L = 200 fim and lj± — 2tt x 1000 Hz. 
The top row shows densities (lighter is higher) of ip+ ( a )j 4>o (b), and ip- (c) for TV = 10 6 atoms, quench time tq = 1 s, and 
final magnetic field strength B/Bq — 1. The vertical lines at ti and ts correspond to the two phase transitions from Fig. [Ha). 
In (d) and (e), faster quenches with tq = 20 ms are shown. The number of atoms in (e) is increased to N — 20 x 10 6 and 
final magnetic field is (B/Bo) 2 = 0.49 to show a cleaner process with a single phase transition and without interaction of fully 
formed domains. Frame (f) shows the number of spin fluctuations in function of time for the cases shown in (d), light line, and 
(e), dark line. 



, 1/3 i 

t ~ Tq , see Fig. E]^a) . However, when the fluctuations 
become sufficiently large, the nonlinearity sets in, which 
is visible as the formation of narrow bubbles of the ipo 
component, see Fig. EJb). This can be seen as the cool- 
ing of the system to the true ground state, which is now 
the phase separated 2C + po state. As the bubbles of 
ipo grow, they narrow to keep the magnetization approx- 
imately conserved in their own neighborhood. We now 
count the number of small bubbles 181, which scales in 



B - Si ~ t/TQ ~ T X 



-2/3 



this fraction is small and scales 

-2/3 



the same way as the number of fluctuations in the Kibble- 

— 1/3 

Zurek theory, Nf ~ Tq .At this point the KZ mecha- 
nism is completed, but it turns out to be just a preludium 
to the ultimate post-selection mechanism that sets the fi- 
nal density of domains Nj. 

After a relatively short time of the growth, as com- 
pared to t, the density of ipo in some of the bubbles 
manages to reach its maximum value no ~ n, forming 
an array of po domains. Their width is equal to several 
spin healing lengths £ s = (B / B) 2 h/ y/2mc2n, which is 
a natural scale defining the minimal width of a po do- 
main, see Fig. SJc). However, it turns out that some of 
the bubbles must have disappeared to keep the magne- 
tization conserved. Indeed, it can be shown that in the 
2C + po ground state the fraction of the system occu- 
pied by the po phase is equal to xo — 1 — B±JB. As we 
are still very close to the critical point, near B such that 



with tq as xq = 1 —B\JB ~ t q 10 - Since the spatial size 
of the bubbles is of the order of several £ s , independently 

2/3 

of tq, wc now have Nd ~ Tq stable domains. This is 
the final post-selected density of domain walls compatible 
with the conserved magnetization. 

Once the domains are stabilized they gradually grow 
in size with xo = 1 — Si/S as the magnetic field is fur- 
ther increased, see Fig. 0Jd). In this last long stage of 
the evolution the domains are stable except for possi- 
ble phase ordering kinetics [l^|. This is a slow process 
where fluctuations eventually merge some domains slowly 
reducing the domain number N^. The fluctuations are 
stronger for a lower number of atoms. In Figure we 
plot the number of created defects for a lower (realis- 
tic) number of atoms N = 10 6 [lj|. The interactions 
are much weaker in this case, which results in a larger 
spin healing length £ s , which sets the smallest possible 
domain size and is responsible for the saturation of the 
number of domains for small tq. Nevertheless, we can 
still clearly see the two different scaling laws for Nf and 
Nd in a wide range of quench times tq. Here, the scal- 
ing exponents appear to be slightly smaller than in the 
previous, "cleaner" case. Since the final state with multi- 
ple domains is not a true ground state, but a metastable 
state |l5} . we attribute this difference to the phase or- 
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FIG. 3. Averaged number of spin domains after the quench 
as a function of the quench time for N = 20 x 10 6 . The scale 
is logarithmic on both axes and Na is decreased by two to ac- 
count for the ground state phase separation into two domains. 
The points are results of truncated Wigner simulations aver- 
aged over 100 runs. The dashed line is the fit to the power law 
with scaling exponent rid = —0.67 ± 0.01. The inset shows 
the maximal number of spin fluctuations counted just after 
the phase transition. The scaling exponent is here found to 
be n f = -0.32 ± 0.01. 
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FIG. 4. Profiles of the density \ip \ 2 for tq = 100 ms, show- 
ing four consecutive phases of the domain formation process 
after crossing the first critical point, (a) At t = 52.5 ms the 
small spin fluctuations begin to grow exponentially, (b) At 
t = 56.25 ms the fluctuations transform into narrow (several 
£ s ) bubbles of ^o- (c) At t = 61.25 ms, as the bubbles mature, 
the post-selection eliminates some of them to keep the mag- 
netization conserved, (d) At t = 100 ms, the domains have 
gradually increased in size and occupy half of the available 
area at B — 2Bi. The formation of the stable po-bubbles 
in the stages (a,b,c) takes place in the narrow time interval 
t = 52.5. ...61.25 around t. The magnetic field B within this 
relatively short time span satisfies B — B\ ~ £/tq ~ rZ 2 ^ . 



FIG. 5. Similar as in Fig. but for a realistic case of N = 10 6 . 
Here, the horizontal dashed line shows the saturation of the 
number of domains due to the finite spin healing length. The 
scaling exponents in this case are rid = —0.71 ± 0.03 and 
n f = -0.35 ±0.01. 

dering kinetics, which effectively decrease the scaling ex- 
ponents [ijj. However, the slow phase ordering taking 
place on a time scale much longer than t should be dis- 
tinguished from the post-selection that happens near the 
same t as the KZ mechanism. 

In conclusion, we investigated a nonequilibrium phase 
transition in a relatively simple and stable experiment 
in an antifcrromagnctic Bosc-Einstcin condensate. We 
demonstrated the occurrence of two scaling laws describ- 
ing the number of spin fluctuations and the final number 
of spin domains. The occurrence of two scaling laws was 
explained in a model of system dynamics including an ef- 
fective post-selection of spin domains due to the conser- 
vation of magnetization. The post-selection transforming 
the scaling law is a general mechanism that should be ef- 
fective whenever the standard Kibble-Zurek mechanism 
is not compatible with an additional conservation law. 
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